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We compare two theoretical approaches to dielectric diblock copolymer melts in an external electric field. 
The first is a relatively simple analytic expansion in the relative copolymer concentration, and includes the full 
electrostatic contribution consistent with that expansion. It is valid close to the order-disorder transition point, 
the weak segregation limit. The second employs self-consistent field (SCF) theory and includes the full electro- 
static contribution to the free energy at any copolymer segregation. It is more accurate but computationally more 
intensive. Motivated by recent experiments, we explore a section of the phase diagram in the three-dimensional 
parameter space of the block architecture, the interaction parameter and the external electric field. The relative 
stability of the lamellar, hexagonal and distorted body-centered-cubic (bcc) phases is compared within the two 
models. As function of an increasing electric field, the distorted bcc region in the phase diagram shrinks and 
disappears above a triple point, at which the lamellar, hexagonal and distorted bcc phases coexist. We examine 
the deformation of the bcc phase under the influence of the external field. While the elongation of the spheres is 
larger in the one-mode expansion than that predicted by the full SCF theory, the general features of the schemes 
are in satisfactory agreement. This indicates the general utility of the simple theory for exploratory calculations. 
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I. INTRODUCTION 

Block copolymers (BCP) consist of several chemically dis- 
tinct sub-chains. They are not only interesting as a model 
system for self-assembly, but also for their chemical versa- 
tility and affordability which have enabled their use in appli- 
cations such as photonic waveguides |1], tough plastics |2], 
ordered arrays of nano-wires |3], etc. At a given chemical 
architecture and temperature, there is one thermodynamically 
stable meso-phase, with typical length-scales comparable to 
the chain size (~ 10-500 nm). However, the material is rarely 
perfectly ordered, but rather is composed of many randomly 
oriented grains of size ~ 1 jim. This has an adverse effect on 
nanotechnological applications. 

A useful way to achieve improved long-range order is to 
subject the BCP sample above its glass transition to an exter- 
nal electric field Erj. Due to the coupling between the field 
and the spatially-varying dielectric constant Jc(r), there is a 
preferred orientation of t he g rains with respect to the field 
HBIlllllJinSIIIlIllIIj- It has been shown by Amund- 
son et al. |4j |5J that the electrostatic free-energy penalty as- 
sociated with dielectric interfaces which are not parallel to 
the electric field direction is the driving force for structures 
to reorient so that their interfaces are parallel to the field 
(Vjc(r) perpendicular to Eq). While the free energy penalty 
can be eliminated by this reorientation of lamellae and cylin- 
ders, it cannot be eliminated in the body-centered-cubic (bcc) 
phase, but only reduced by distorting the bcc spheres. Thus, 
the free energy of this distorted bcc phase, whose symmetry 
is reduced to R3m, increases with respect to the full disor- 



dered liquid (dis), lamellar (lam) and hexagonal (hex) phases 
lUlll . a circumstance which can bring about a phase transition. 
The effect of the electric field on the BCP morphology has 
been substantially accounted for recently 11411 by incorporat- 
ing the electrostatic Maxwell equations in the full set of self- 
consistent field (SCF) equations, which permits calculation of 
the phase diagram at arbitrary degrees of segregation. 

In this paper we compare two theoretical approaches to 
such a system; the aforementioned SCF study and a simple 
analytical approximation consisting of a Ginzburg-Landau ex- 
pansion of the free-energy 1 15], valid only close to the order- 
disorder temperature (ODT). The paper is organized as fol- 
lows. In Section 2 we present the free-energy model which 
includes the electrostatic energy of the BCP in the field. In 
Section 3 we calculate the way in which an initial mesophile 
deforms under the influence of the field, and also find the rel- 
ative stability of the competing phases. A comparison is made 
with the results of the SCF model. Section 4 contains a brief 
conclusion. 



II. MODEL 

Although the effect we consider here is generic to any 
multi-block BCP melts, we will restrict the discussion in this 
paper to the simplest A/B af/-block copolymer, where a spatial 
variation of the relative A/B monomer concentration yields a 
spatial dependence of the dielectric constant and, hence, of the 
response to an external electric field. 

We also assume for simplicity that the A monomeric vol- 
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ume is equal to the B one. Then the volume fraction of the 
A monomers /, (0 < / < 1), is equal to its molar fraction. 
The order parameter <p (r) is defined as the local deviation of 
the A-monomer concentration 0a (r) from its average value: 
<p (r) = <pA (r) — /■ From an incompressibility condition of the 
melt we also have at each point r, <p B (r) = 1 — 0a (r). In the ab- 
sence of any external electric fields, the bulk BCP free-energy 
per polymer chain, F b , in units of k B T, can be written as a 
functional of the order parameter, <p(r). One way of generat- 
ing a simple analytical expansion in the order parameter relies 
on a Ginzburg-Landau-like free energy, which can be justified 
close to the order-disorder point (ODT) Il6l fl7ll and is 
repeated here without further justifications: 
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b is the Kuhn length, R g the radius of gyration, % is the Flory 
parameter, N = Na + Nb the total chain polymerization index, 
NXs the spinodal value of j/V 1 15], c is a constant of order 1, 
and A and u are functions of / as in refs. 0E0. The 
phase diagram in the ( f,%N) plane, as derived from the free- 
energy, eq[2 is symmetric with respect to exchange of / and 
1 — /. For small values of % ~ l/T, the melt is disordered: 
(r) = / is constant. For %N larger than the ODT value of 
~ 10.5 and for nearly symmetric BCP composition (f ps j), 
the lamellar phase is the most stable. As \f — i | increases, 
the stable phases are doubly-connected gyroid, hexagonal and 
bcc phases 1 1 5l M 61. 1 1 t SII . 

Let us now consider a BCP slab placed in an external elec- 
tric field, E<). The free-energy per polymer chain, again in 
units of k B T, is F m = Ft,+F &s , where the electrostatic energy 
contribution F es is given by the integral 
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Here Eq is the vacuum permittivity, >c(r) is the local dielectric 
constant, v p is the volume per chain, and y/ is the electro- 
static potential obeying the proper boundary conditions on the 
electrodes. The local field is E(r) = — Vy/. We note that the 
variation of F es with respect to y/ yields 



V-(eofc(r)Vv/)=0 



(4) 



BCP slab filling the gap between two parallel and flat elec- 
trodes separated by a distance d and potential difference V. 
Even when a non-homogeneous dielectric material like a BCP 
fills the gap between the two electrodes, the spatially averaged 
electric field in between the electrodes (E) is constrained to be 
Eq = V jd. The local field E(r) differs from its average due to 
the nonuniformity of the dielectric constant, since K = ff(0) 
depends on the local concentration 0(r) through a constitu- 
tive equation. In this paper we assume for simplicity a linear 
constitutive relation, 



jc(r) = (jc) + 0(r)Ajc, Ak=k a -k b , 

(k) = fK A + (l-f)K B , 



(5) 



where throughout this paper we use K A =6.0 and K B 
2.5, thus modelling an A/B diblock copolymer where the A 
block is polymethylmethacrylate (PMMA) and the B block is 
polystyrene (PS), as is used in many experiments. Other con- 
stitutive relations can be considered 1 19]. 

When a field is applied on a melt in the lamellar or hexago- 
nal phases, it exerts torque which causes sample rotation. The 
torque is zero, and the energy lowest, when the lamellae or 
cylinders are oriented parallel to the field. In such states, as 
well as the disordered phase, the electrostatic energy, eq[3] of 
the system is equal to a reference energy, given in eq[6] The 
bcc array of spheres, on the other hand, always has dielectric 
interfaces that are not parallel to the field, and its electrostatic 
energy is higher than the reference value. Hence, the spheres 
elongate in the applied field direction, to an extent which is a 
balance between electrostatic and elastic forces, as calculated 
below. 

The reference energy per polymer chain, in units of k B T, is 
simply — j(k)Eq, where Eq is the applied field measured in 
the natural unit (k B T / Eqv p ) 1 I 2 , 
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Let us estimate the value of the actual applied field cor- 
responding to Eq = 1. At 100° C and using typical polymer 
volume per chain in the range v p ~ 50 — 250 nm 3 , we find 
Eq ~ 47 — 107 V/jim. This is a relatively large field that can 
cause dielectric breakdown in some BCP films. Therefore, the 
experimentally interesting regime is usually Eq < 1 . 

The free-energy F tot as formulated above is valid close to 
the ODT point (weak segregation limit), where the concentra- 
tion variations are small, 0(r) 1, and therefore the analysis 
can be carried out within the so-called one-mode approxima- 
tion. Motivated by recent experiments 00, we concen- 
trate on the transition from distorted spheres to cylinders or 
disordered melt in presence of an applied electric field. Taking 
Eo to be in the (1,1,1) direction, we write the order parameter 
<p as a linear superposition of six components 
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where 



which is the usual Maxwell equation V-D-0 for the displace- 
ment field D = EqkE. We consider a simple geometry of a 
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The q's and k's are wave-vectors given by 



and 
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and all have the same magnitude go- The three linearly depen- 
dent q, are orthogonal to the (1,1,1) direction and describe a 
hexagonal phase with axis along that direction. The three k, 
have equal and non-zero projections on the (1,1,1) axis. The 
six wavevectors transform into one another under the symme- 
try operations of the bcc phase. In the absence of an exter- 
nal field, each of these wavevectors would contribute equally 
in the order parameter expansion, so that g and w would be 
equal. These wavevectors characterize the first mode in such 
an expansion. Hence the name of the approximation. 

The amplitudes w(Eq) and g(Eo) depend on the magnitude 
of the average external field Eq. Depending on the values of 
the two amplitudes, g and w, we can represent the order pa- 
rameter of all phases of interest in the form of eq0 w = ^0 
represents an undistorted bcc, while an R3m (distorted bcc) 
phase oriented along the (1,1,1) direction is represented by 
two non-zero amplitudes w ^ g. A hexagonal phase of cylin- 
ders whose long axis is in the (1,1,1) direction has w^O and 
g = 0. And finally, g = w = represents the disordered melt. 
As was mentioned above, the spatially averaged electric field 
is simply the magnitude of the external field, Eq. However, 
local changes in 0(r) give rise to local changes in K"(r). Con- 
sequently, the electric field can be written as follows: 



E = E + c5E 



(10) 



where 8E is the deviation from the average. The symmetry of 
the electrostatic potential follows from the Maxwell equation, 
eq. 0] the constitutive equation, eq. |5J relating the dielectric 
constant to the order parameter, and the R3m symmetry of the 
distorted bcc phase. From the symmetry of the potential, one 
finds that 8E can be decomposed in terms of the same bases 
of k's and q's given above: 



8E = aEi+j3E 2 

3 

Ei = £b£q;cos(q r r), 

3 

E 2 = E £k\cos(k r r) (11) 

i=i 

where k; = k;/|kj| and q,- = qj/|q/| are unit vectors. Within 
the one-mode approximation, the local field at each point is 
given in terms of two amplitudes, a and j3 . They determine 
the projection of the deviation field SE onto the Ei or E 2 di- 
rections. 

In order to proceed we insert the E'-field expressions, ea ll II 
into the electrostatic free-energy, eq[3] Using the definitions 
of eqsEEJand the properties: 

<ln ; qrq; = -<ZoA ' •' ./ 
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(1,1,1) ; krq^O ; k i -q j = ±^, i^j (13) 



we can perform the rather straightforward spatial averages of 
the various terms in the free energy, eq|3] and obtain: 
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The last term is simply the reference energy which is common 
to all phases. 

For a given state of (a given BCP morphology), which 

is determined by a given value of w and g, the values of 

a(w,g) and P(w,g) are determined by the Maxwell equation, 

eq |4] This is equivalent to obtaining them by varying F ss 

with respect to a and j3. The procedure yields a = and 
l 



j8 = - y/2/3AKg / ((k) + jAkw). Thus, F es is given by 



(Ak) 2 



2(k) + wAk 



2 X ' 



(15) 



It is instructive to compare this result with the perturbation 
expression used by Amundson, Helfand and co-workers | 
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a result obtained from eq^]in the limit Ak/(k) — > 0. In 
particular, the result of eq ^j] yields a free energy which is 
symmetric under the interchange of monomers A and B, while 
the result of eq|15ldoes not. Such a symmetry is not expected 
in general: a system of ultra-high K spheres (metallic limit) 
embedded in an insulator matrix has a different energy than 
the system of insulating spheres embedded in a ultra-high k 
matrix, even when the average dielectric constant is the same 

(K). 

Employing the single-mode Ansatz = w0i +g<p2 in e q[D 
we finally obtain for the total free-energy per polymer chain 
in units of kgT 
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In the next section we minimize this energy with respect to 
w and g at a given dimensionless external field Eq and polymer 
architecture /, calculate the elongation the spheres of the bcc 
phase, and obtain the phase diagram. 
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FIG. 1: Normalized free-energy per polymer chain /„, defined in eq 
1181 of the distorted bcc phase (R3m) as a function of dimensionless 
field Eq. The system is characterized by / = 0.3 and %N = 14.4. 
We compare the one-mode calculation (solid line) as obtained from 
minimization of eq |17l with a SCF calculation (dashed line). The R3m 
phase in the SCF calculation has a lower free energy that the solid 
line (one mode), and crosses the hex energy at higher value of Eq of 
about 0.49, while the one-mode approximation crosses at Eq ~ 0.43 
(both marked with arrows). In this figure and following ones we used 
K A = 6 and Kb = 2.5, modelling a PMMA-PS copolymer. 
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FIG. 2: The amplitudes w and g normalized by their common value 
at zero Zs-field, g(0) = w(0), as a function of dimensionless exter- 
nal field Eq. Solid line: one-mode approximation. The amplitudes 
have a discontinuous jump at Eq — 0.43, where the structure contains 
cylinders oriented along the field (g = 0). Dashed lines: the same, 
but taken from a multi-mode SCF calculation with a jump at higher 
E values of Eq ~ 0.49. All parameters as in figure 1. 



III. RESULTS 

As noted above, the functional form (j) = wty\ + g(j>2 allows 
us to distort smoothly a bcc array of spheres (having nonzero 
w — g) via a distorted bcc phase (w ^ g), and into a hexagonal 
array of cylinders (non-zero w but with g = 0). The disordered 
phase is given by w = g = 0. One is thus able to obtain the 
full phase diagram by minimizing eq with respect to the 
amplitudes w and g. 

Before presenting the phase diagram, let us consider a point 
in the (f,Nx) plane for which the stable phase at zero E field 
has a bcc symmetry. For presentation purposes, in figure 1 
we have subtracted from the free energy the reference elec- 
trostatic energy, — (k)Eq/2, common to all phases, also sub- 
tracted the total free energy of the bcc phase in zero field, 
f t Q( C (0), and normalized the resulting free energy by that of 
the hex phase in zero field; that is we have plotted 

_ F tot (E (i ) + (K)E^/2~F^{0) 
fn ~ (0) 

In the figure we show how the free energy /„ changes with 
E for / = 0.3 and N% = 14.4. At Eq = the bcc is the sta- 
ble phase, and its free energy increases with increasing field 
Eq, until it equals the free energy of the hex phase at a tran- 
sition field Eq ~ 0.43. At larger fields the stable structure 
is a hex phase of cylinders oriented along the external field 
Eq. The solid line in figure 1 is the result obtained from the 
one-mode approximation given above, while the dashed line 
is a obtained from the SCF theory, (as in Ref ffl v It has 



a lower free-energy. Consequently, the transition field in the 
SCF framework is higher and occurs at about Eq ~ 0.49. 

Figure 2 is a plot of the amplitudes w(Eq) and g(Eo), nor- 
malized by their zero-field value w(Eq = 0) = g(EQ = 0). Both 
amplitudes start at their common value in the undistorted bcc 
phase. As the field increases, w increases while g decreases. 
The spheres elongate in the direction of the field as a result 
of competition between electrostatic and elastic forces. At the 
transition field, there is a sharp, discontinuous, transition in 
the order parameter. Above this field, w attains a fixed value 
while g drops abruptly to zero. In this state the BCP mor- 
phology is that of cylinders oriented parallel to the external 
field. The dashed lines correspond to the values obtained from 
the SCF theory. Clearly, in the one-mode approximation, the 
spheres' deformation and eccentricity are larger than in the 
SCF theory. 

The above calculation can be repeated for any {f,N%) 
and Eq field values and allows the construction of the full 
three-dimensional phase diagram in the (/,Nx,Eq) parame- 
ter space. In figure 3 we present a cut of the phase diagram 
at fixed / = 0.3. The region of a stable R3m phase (distorted 
bcc) is bound by two lines of phase transitions: one between 
this phase and the disordered phase, and the other between it 
and the hex phase. These two lines meet at the triple point 
(Xt,E t ). In figure 3, the different triple point values obtained 
from the two calculations are used to rescale both axes: x/Xt 
and EQ/E t . At fields larger than E t the R3m is not stable at 
any value of X- The solid lines in the figure are the one-mode 
prediction, while the dashed lines are obtained with the SCF 
calculation. The values of E t are 0.49 and 0.67 for the two 
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FIG. 3: Phase diagram in the (x/%t< *o/Et) rescaled plane for a fixed 
value of BCP asymmetry, / = 0.3. The distorted bcc (R3m) region 
is confined between two transition lines which terminate at a triple 
point (Xi,E r ). The other two phases are the disordered one (dis) and 
the hexagonal one (hex). The solid lines are the prediction of the 
one-mode approximation, with axes scaled by the appropriate Xt — 
143 /N and E t ~ 0.49. Dashed lines are the SCF calculation scaled 
by the appropriate SCF values of % t ~ 14. 1 1 /N and E, ~ 0.668. 



theories, respectively. 

For an additional comparison between the theories, we have 
examined the case in which, at a fixed value of / = 0.3, the 
dielectric constants of the majority and minority components 
are interchanged (i.e. = 6.0 <-*■ Kb = 2.5, hence, AfC — > 
—Ak). In both theories we find an increase in the value of 
the external field needed to bring about a transition from the 
distorted bcc phase to the hex phase. Thus, this subtle effect, 
which is not captured by the perturbation result of eq^] is 
obtained in the simple one-mode approximation, ea!15l 

IV. CONCLUSIONS 

A simple theory for a non-homogeneous diblock copoly- 
mer (BCP) melt in an external electric field is presented, and 
compared with a more accurate, but more computationally in- 
tensive self-consistent field (SCF) one. The differences be- 



tween the two theories in zero external field are well known. 
In particular, the accuracy of the phase boundaries produced 
by the one-mode approximation deteriorates outside the vicin- 
ity of the ODT point (weak segregation) as compared to the 
SCF theory 1 18]. However, as in the zero field case, the qual- 
itative behavior of the system in the presence of a field is de- 
scribed surprisingly well. The full electrostatic free-energy 
contribution is used, consistent within the one-mode approxi- 
mation, eq^l This was not accounted for in previous analyt- 
ical studies HIM ITU Il3ll . where only quadratic terms in the 
electrostatic potential were retained. The simple one-mode 
approximation captures the elongation of the spheres of a bcc 
phase when placed under an external field. The elongation is 
in the direction of the applied £0 field. The two amplitudes 
describing this elongation, w and g, as shown in figure 2. At a 
threshold value of the electric field, a first-order transition to 
a hexagonal phase occurs and the amplitudes jump discontin- 
uously. 

As shown in figure 3, the simple, analytic, one-mode ap- 
proximation also captures the essence of the phase diagram; 
the reduction in the phase space occupied by the distorted bcc 
phase as the field increases, and its eventual disappearance at 
a triple point. 

Lastly the one-mode theory also captures the subtle inter- 
play between structure and electrostatic response as evidenced 
by its prediction of a different critical field for phase transi- 
tions when the dielectric constants of the constituents are in- 
terchanged, a prediction in agreement with the more accurate 
theory 0. 

Given its ability to capture all of the above effects, and 
given its extreme simplicity, such a theory could serve for use- 
ful exploratory studies in other problems concerning the effect 
of electric fields on block copolymers. 
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